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We calculate the quark number susceptibility in the deconfined phase of QCD using the hard ther- 
mal loop (HTL) approximation for the quark propagator and quark-meson vertices. This improved 
perturbation theory takes into account important medium effects such as thermal quark masses and 
Landau damping in the quark-gluon plasma. We explicitly show that the Landau damping part 
in the quark propagator for space-like quark momenta does not contribute to the quark number 
■ susceptibility due to the quark number conservation. We find that the quark number susceptibility 

' only due to the collective quark modes deviates from the free one around the critical temperature 

but approaches free results at infinite temperature limit. The results are in conformity with recent 
lattice calculations. 
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In recent years substantial experimental and theoretical efforts have been undertaken to investigate the versatile 
physics issues involved in ultra- relativistic heavy-ion collisions, i.e., collisions of atomic nuclei in which the center- 
of-mass energy per nucleon is much larger than the nucleon rest mass. The principal goal of this initiative is to 
I explore the phase structure of the underlying theory of strong interactions - Quantum Chromodynamics (QCD) - by 
psj ■ creating in the laboratory a new state of matter, the so-called Quark-Gluon Plasma (QGP). This new state of matter 
' is predicted to exist under extreme conditions like at high temperatures and/or densities, when a phase transition 
takes place from a hadronic to a deconfined state of quarks and gluons Such information has essentially been 
confirmed by numerical lattice QCD calculations at finite temperature, which show a rapid increase in energy 
density and entropy density as a function of temperature. Numerical solutions of QCD also suggest that the critical 
' temperature is about 170 MeV and provide information on the equation of state [Q. 
i-S^ ' The various measurements taken at CERN SPS within the Lead Beam Programme do lead to strong 'circumstantial 
^pH, evidence' for the formation of the QGP (|J^. Evidence is circumstantial as any direct formation of the QGP cannot 
O ^ be identified. Only by some noble indirect diagnostic probes like the suppression of the J /'^ particle, the enhanced 
D , production of strange particles, especially strange antibaryons, excess production of photons and dileptons, the forma- 
' tion of disoriented chiral condensates, etc. the discovery can be achieved. An extensive amount of theoretical study 
has also been devoted over the last two decades in favour of these well accepted probes of a deconfincment (QGP) 
phase. 

Recently screening and fluctuation of conserved quantities have been considered as important and relevant probes 
^ , of the QGP formation in heavy- ion collisions 0-0- In the confined/chirally broken phase charges are associated with 
the hadrons in integer units whereas in the deconfined/chirally restored phase they are associated with the quarks in 
fractional units which could lead to charge fluctuations which are different in the two phases ■ The fluctuations 
can generally be related to the associated susceptibilities [0,0- The quark number susceptibility is associated with 
the number fluctuation which measures the response of the number density with infinitesimal change of the quark 
chemical potential. Hence the quark number susceptibility can be related to charge fluctuations and is therefore of 
direct experimental relevance. The quark number susceptibility has been investigated in lattice QCD simulation |l^ ] 
which showed that it is zero at low temperature and rises suddenly to nonzero values across the deconfinement phase 
transition. At high temperature QCD it has been analysed [|l3[ to show non-perturbative temperature effects at 
next-to-leading order. Recently, it has been discussed |Q in connection with the role of the fluctuations during 
the dense stages of the collision to exploit the electromagnetic probes with the hadronic probes. A very recent 
lattice simulation |T5| has verifled a new relation between susceptibilities and screening masses and explains that the 
non-perturbative phenomena are closely connected with deviations from weak coupling limit or bare perturbation 
theory, indicating the need to resum the weak coupling series. The purpose of the present calculation is to investigate 
the quark number susceptibility within the HTL-resummed perturbation theory which incorporates non-perturbative 
effects such as effective masses of the collective quark modes (quark and plasmino modes in the medium originating 
from the poles of the HTL-propagator) and Landau damping for space-like quark momenta, reflecting the physical 
picture of the QGP as a gas of quasiparticles. As we will see below, the quark number susceptibility obtained in HTL 
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approximation is in agreement with recent lattice [ p^ observations. 



A. Fluctuation and Susceptibility: 

Let Oa be a Heisenberg operator. In a static and uniform external field J-'a, the (induced) expectation value of the 
operator Oa (0, ~af ) is written as 

Tr to fO ^U~'^('^+'^'=="'] 1 r 
K - {Oa (0,1?))^ = ^ Tr ^ ^ F 7 ^"'"^^^ ' 

where translational invariance is assumed and Tiex is given by 

ne. = -J2 f d?x Oa (0, a*) . (2) 

The (static) susceptibility Xap is defined as 

Xaa{T) 



p [ Sx (O„(0,x^)O, (0, O')) , (3) 
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assuming no broken symmetry (Oq (0, Ir^)) = (Oo- (O, O)) =0. (0^(0, x)C'ct(0, 0)) is the two point correlation 
function with operators evaluated at equal times. 



B. Quark Number Susceptibility: 



The quark number susceptibility is the measure of the response of the quark number density with infinitesimal 
changes in the quark chemical potential jiq + 5^,q. Under such a situation the external field, Ta, in (^ can be 
identified as the change in quark chemical potential and the operator Oa as jo = 97o9i where jfj_{t,x) — qj^q is 
the vector meson current. Then the quark number susceptibility for a given quark fiavour follows from (0) as 
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13 j d'x (jo(0,x)jo(0,0)) =/3 / d^x Soo{0,x-) 
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where 5*00(0, x) is the time-time component of the vector meson correlator S^^{t, x) — {jfj.{t, x)ji,{0, 0)) and the number 
density can be written as 



Pq 



1 Tr [Afq, 



V Tr re-/5C«~A',A/',)l 
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1 dn 
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with the quark number operator, J^q ~ J jo{t,x) d^x, and fl = — TlnZ is the thermodynamic potential and Z the 
partition function of a quark-antiquark gas. 

Taking the Fourier transform of 6*00(0, x), it can be shown that 



Xq (T) = lim P 

p— »o 



diu 
2^ 



^oo i(^,p) 



Using the fluctuation-dissipation theorem it can further be shown that [|ri[ 
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^ /•+°° dcj 



■Imnoo (cj,0) , 



where n^.(w,p) = FT -i0it) ( j^(<, x), jV(0, 0) 



and FT stands for Fourier transformation. 
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1. Free case 
K 




Q = P-K 



Fig. 1 The self-energy diagram for free quarks. 

To lowest order in perturbation theory one has to evaluate the time-time component of the self energy diagram 
shown in Fig 1, where the internal quark lines represent a bare quark propagator Sf{L) which can be expressed in 
the helicity representation for massless case as [L — {la, I), I — [0 

with 

d±{la,l)^-h±l . (9) 
The corresponding spectral function is given by 

p£(/o,0 = '5(?otO- (10) 

Now the time-time component of the vector meson self-energy in Fig 1 can be written as 



ko ^ ' 



(11) 



where, Q = P — K, and Nf and Nc are, respectively, the number of quark flavours and colours. 

Substituting the propagator (||) in (pi]), performing the trace, and following the summation relation of Rcf. [ p7[ , 
we extract the imaginary part as 

lmn™(a;,p = 0) = 4iV/7Ve7r (l - e'^") J -0^ J dx J dx' 5 {u - x ^ x') 

np (x) np {x') pi^ (x, k) pL {x , k) , (12) 

where x and x' are the energies of the internal quarks and np is the Fermi distribution function. Using ( p^ in ([T^), 
we get 

/d'^k 
^np (k) {1 - np (k)) . (13) 
{2ir) 

Inserting ( p^ in (|^), we obtain the quark number susceptibility in lowest order perturbation theory |^,|ll|] 

/d^k e'"^ 

Alternatively, considering the lowest order thermodynamic potential of a quark-antiquark gas |Iq,M] 
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~2NfNcT 



(2^ 



/3Ek + In (l + g-ZSt^^^-A',)) + In (i + e-P^^>'+^^^^ 



(15) 



one could also arrive at the same expression using (Q) and (g). 

Due to the quark number conservation Im n(a;, 0), as obtained in (p^, is proportional to S{llj). This leads to a 
general relation between the quark number susceptibility and the time-time component of the polarisation tensor in 
the vector channel 



xl (T) = -mfN^I ^ ^ noo(o,o) ^ 



{2T:y dk 



(16) 



which provides a connection | pd] , p^ , p^ -|20|] to the electric screening mass, This relation is only valid in lowest 
order in perturbation theory. 



2. HTL case 



Now we turn to the estimate of the quark number susceptibility beyond the free quark approximation by invoking 
the in-medium properties of quarks in a QGP. In the weak coupling limit (5 <C 1), a consistent method is to use 
the HTL-resummed quark propagators and HTL quark- meson vertex if the quark momentum is soft (~ gT). Using 
this improved perturbation theory at least to some extent non-perturbative features of the QGP such as effective 
quark masses and Landau damping are incorporated through the effective quantities like quark propagators and the 
quark-meson vertex. 

Let us, however, note that we do not aim at a complete leading order perturbative calculation. Rather we want 
to study the influence of medium effects incorporated in the HTL resummed quark propagator. Hence, we will use 
this propagator for the entire momentum range instead of consistently distinguishing between soft and hard momenta 
pl| . Anyway, since this distinction is only possible in the weak coupling limit, (7 <C 1, it cannot be used in our case, 
in which we want to compare our results to QCD lattice calculations. The approach, considered here, is in the same 
spirit as the one for calculating meson correlators in the QGP It should be noticed that our results are gauge 
independent due to the gauge invariance of the HTL quark propagator. 

The HTL-resummed quark propagator, S*{L), can be obtained [0| from (@) by replacing d±{L) as 



D± ilo,l) = -lo±l 
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Qo 



TQi 



(17) 



where the thermal quark mass is given by — g(T)T/\/6, and Qn{y) is the Legendre function of second kind. The 
HTL- vertex can be obtained lESl as 



(18) 



with 



dn 



All {R-Pi){R-P2) 



^G^^ (-Fi,-F2) , 



(19) 



where R = (— l,r) is a light-like four vector, R^ = 0. The effective propagator and vertex are related via the Ward 
identity. 

The HTL-spectral function reads jl^ 



p±(/o,0 = 



;2 _ /2 
2m^ 



[S{lo - uj±) + S{lo + u;^)] + f3±{lo,l)e{l^ - ll) 



(20) 



with 



±lo-l 



l{-lQ±l)+ml (±1 



21 l-lo 



"q I 



(21) 



Here the zeros uj± (l) of D±{L) describe the two branches of the dispersion relation of collective quark modes in a 
thermal medium [|T^ . Furthermore the HTL-resummed quark propagator acquires a cut contribution below the light 
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cone {Iq < as the quark self-energy has a non- vanishing imaginary part, which can be related to Landau damping 
for space-like quark momenta resulting from interactions of valence quarks with gluons in the thermal medium. In 
addition, an explicit temperature dependence only enters through mq(T) as well as through the strong coupling 
constant, g'^{T) = 47ras(T) with 

12tt 

"^^^^ ^ (33-2iV/)ln(QVA2)' ^^^^ 

where Aq = 200 - 300 MeV. For the momentum scale Q we take the energy of the lowest Matsubara mode Q — 2ttT . 
For checking the sensitivity of the susceptibility to uncertainties in the coupling constant we will use also Q = AttT. 
It should also be noted that the HTL-resummed propagator is chirally symmetric in spite of the appearance of an 
effective quark mass jlTt . 

Now we need to calculate the imaginary part of the time-time component of the self energy diagrams given in 
Fig. 2, in which blobs represent the effective quantities. The tadpole diagram in Fig. 2b is essential to satisfy the 
transversality condition, P^n^''(P) — 0. As will be seen below this has a very important effect on the quark number 
susceptibility by partially compensating the cut contribution in diagram Fig. 2a. 

K 




(a) - - - ^ - 

P \ IP 



Q^P~K 



(b) 



Fig. 2 The self-energy (a) and tadpole (b) diagrams for quarks in the HTL-approximation. 

Now the time-time component of the self-energy in diagram Fig. 2a can be written as 

n?" (P) NfN^TY, J ^Tr [S* {K) pO (K - P, -K; P) S* (Q) (P - K, K; -P)] . 



(23) 



k 

The time component of the HTL- vertex can be obtained by carrying out the angular integration in (|l8|) for ~jf = 

•2 



where 



r'{K~P,K;P)= h- -^SQo j 7" + ^^Q^''-^ ' (24) 



SQn=Qn[^)-Qa(^^) ■ (25) 



k ""V k 



Alternatively, r°(Pi,P2) can also be obtained from the Ward identity P^r''(Pi, Pa; -P) = 5'*"^(Pi) - S'*"^(P2). 
Now performing the traces in (|23|), we get 



nr(Po,#=o)^2iv,iv,r^ / ^ 



(a + 



(a - 6)2 



i?+(if)i?_(g) (i^) (Q) 



(26) 



where 



1 - ^^Qo 
Pofc 



pok 



-^)-«»(f)(-f) 



(27) 



where go — Pa~ and Q„ (— 2/) = (— l)""*^^ (y) has been used in (p7|). Following the summation formula of Ref. |Tl 
the imaginary part of ( p6[ ) then can be written as 

lmn?"(cj,p = 0) = TV/TVeTT (1 - e'^'^) J y" f^a; y (cj - x - a;') «f (a;) (x') 

P+ {x, k) p_ (x', k) 

{i(i-f).,(.',*0 + 1(1 + 1). 
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x + a;' 



4::;^e (fc2 - 



The contribution of the time-time component of the tadpole diagram in Fig. 2b can be written as 
n°°(P) =NjN,Tj2 J j^TT[S*iK)r°''i-K,K;^P,P)] , 

ko 

where the effective HTL four-pomt function can be obtained from the relation 

P^r^^(-X, K; -P, P) = r^(X - P, -K; P) - V'i-K - P, K- P). 
At ^ = the four point function is obtained as 

TTl Tf~l 

r°° i-K, K; -P, P) - - ^ {SQo + 5Q'o) 7° + ^ {SQi +SQ\)k-j, 



with 



SQ'n = Qn 



ko +po 



Proceeding exactly the same way as before, we get the imaginary part of the tadpole 

Imllf{uj,p = 0) = NfNai: (l - e^'^) [ -^-^ [ dx [ dx'S {lo - x - x') np (x) np {x') 

J (27r) J J 



4^9 (fc2 - x') 



1 



- 1-- U+(:r',fc) + - 1 + - p_ (x',fc) 



(28) 

(29) 

(30) 
(31) 

(32) 

(33) 
(34) 



It can be seen that the tadpole contribution compensates the second term of ( [28D and the total contribution becomes 

d^k 



lmn°"(w,p==0) =ANfN^Ti{l-e'^'^) J -^-^ J dx J dx' S {to - x - x') np {x) np (x') 



X + x' 



(2^)3 

I 

p+ (x. k) p- {x\ k) 



(35) 
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Now combining (|2^), ( |35| ) and the quark number susceptibility in HTL-approximation is obtained as 



_fk_ 

(2^ 



K(fc)-fc^)' 

4m^ 



(36) 



Since the quark number susceptibihty is constructed from two quark propagators, in general they should receive 
pole-pole, pole-cut and cut-cut contributions [ pT[ . However, according to ( p6| ) it has only pole-pole contributions from 
the two collective quark modes (quark and plasmino modes in the medium). The cut contributions (pole-cut and 
cut-cut) due to space-like quark momenta (Landau damping) do not contribute because of the number conservation 
in the system. In the high temperature limit, the second term due to the plasmino mode decouples from the medium 
whereas the first term reduces to the free susceptibility given in (p^). 
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Fig 3. Xq{T)/NfT'^ as a function of T/Tc in the HTL-approximation. Nf is the number of quark flavours and Tc is 
the critical temperature for the deconfined phase transition. For the solid curves ( p^ ) with Tc/Aq — 0.49 and for 
the dashed curves Aq = 200 MeV and Tc = 170 MeV have been used. The lower curves correspond to Q = ^ttT and 
the upper ones to Q = AttT. 



We now discuss our results for the quark number susceptibility. In Fig. 3 we present the temperature dependence 
of the quark number susceptibility Xgi'T) '^^ units of NfT'^, the free field theory susceptibility. We have used the 
temperature dependence of the strong coupling constant as given in (|2^) for two different ratios of Tr, to Aq and two 
different values for Q. The choice Tc/Ao = 0.49 has been used in the lattice calculations of Ref. jl^. We observe 
that the quark number susceptibility does not depend strongly on the choice of the coupling constant. (For Tc = 170 
MeV and Aq = 300 MeV the susceptibility decreases by about 5% compared to the dashed curves.) The susceptibility 
increases with the increase of temperature. It is interesting to note that the result lies significantly below the free 
result even at T = 5Tc and the deviation is 25 — 40%. This deviation agrees with the recent lattice observation |l^] 
along with a slow approach to the free quark susceptibility. This result is reminiscent of the fact that the quark 
number susceptibility in HTL-approximation contains non-perturbative information about the QGP phase at high 
temperature. In the infinite temperature limit, the HTL result exactly reproduces the free quark susceptibility. This 
can clearly be seen in Fig. 4 below. 

In Fig. 4 we display the quark number susceptibility in units of NfT^ as function of rriq/T using ( p2| ) with Aq = 300 
MeV and Q = 2ttT. As expected for small thermal quark masses, the HTL-susceptibility begins with the free field 
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theory result, NfT"^, with a flat tangent at niq = and a quadratic dependence on rUq. With the increase of the 
mass the susceptibility starts decreasing monotonically. For low values of the thermal quark mass, the susceptibility 
is large due to the fact that it is relatively easy to create an additional quark or antiquark. On the other hand, if 
quarks acquire thermal masses in the medium, the susceptibility decreases, which could qualitatively be understood 
as an effect of the Boltzmann factor. 
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Fig 4. Xq{T)/NfT'^ as a function of niq/T for free quarks (dashed line) and in HTL-approximation (solid line). rUq 
is the thermal quark mass. 



C. Conclusion: 



We have calculated the quark number susceptibility in the HTL-approximation which incorporates in-medium effects 
of the QGP phase such as quark masses and Landau damping through the HTL resummed propagators and HTL 
quark-meson vertex in the vector meson channel. We have discussed the influence of the various non-perturbative 
effects on the quark number susceptibility in the deconflned phase. The Landau damping contribution due to space- 
like quark momenta drops out because of the quark number conservation in the system. Technically this occurs partly 
due to a cancellation of the two diagrams in Fig. 2 containing HTL-resummed quark-meson vertices, and partly due 
to kinematical reasons in (|35|). 

We find that the quark number susceptibility is significantly smaller than the free susceptibility for moderately high 
temperatures. Our results are in good agreement with recent lattice calculations. Besides the on-set of confinement 
and chiral symmetry breaking close to Tc, HTL effects may explain the lattice results similar as in the case of the free 
energy p6| but in contrast to the meson correlation functions 
to charge fluctuations, it could be an interesting observable in heavy-ion collisions. 

Note added: After completing our investigation a preprint by Blaizot, lancu, and Rebhan appeared on the same topic 
. They computed the quark number susceptibility from the thermodynamic potential, obtained in an approximately 
self-consistent resummation of HTLs. However, they do not use the effective HTL resummed vertices, which plays 
an important role by partly compensating the Landau damping contribution in the quark number susceptibility as 
shown in our investigation. As in the present paper, they found results similar to lattice data. 
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22 1 . Since the quark number susceptibility is related 
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